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Abstract 

In this paper we analyze the rural-urban migration phenomena as it is usu- 
ally observed in economies which are in the early stages of industrialization. The 
analysis is conducted by means of a statistical mechanics approach which builds 
a computational agent-based model. Agents are placed on a lattice and the con- 
nections among them are described via an Ising like model. Simulations on this 
computational model show some emergent properties that are common in develop- 
ing economies, such as a transitional dynamics characterized by continuous growth 
of urban population, followed by the equalization of expected wages between rural 
and urban sectors (Harris- Todaro equilibrium condition), urban concentration and 
increasing of per capita income. 

1 Introduction 

Economic development generates significant structural transformations, such as changes 
in the demographic condition and the production structure. As mentioned in chapter 10 
of ref. P] the most important structural feature of developing economies is the distinction 
between rural and urban sectors. The agriculture plays a key role in the development of 
the urban sector of an economy: it provides both the food surplus that enables the urban 
sector to survive, and the supply of labor to the expanding industrial urban sector. As 
suggested in chapter 5 of ref. , the fundamental part of the transformation from mostly 
dispersed rural and agrarian countries in a more concentrated urban and industrial is the 
flux of a large number of individuals through migration from rural areas to urban areas. 
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In this paper we examine the rural-urban migration phenomena which takes place 
during the industrialization process. The analysis is carried out by using an agent-based 
computational model which aims to describe some of the main structural features of a 
developing economy. We look at the rural-urban migration as a discrete choice problem, 
which allows us to formalize the migration process by using an Ising like model. We 
modelled the migratory decision in the usual manner considering the force exerted by the 
difference of earnings between the sectors. Moreover, it is also included a new factor, the 
influence from neighbors like in the Ising model. 

The paper is organized as follows. Section |21 describes the economic setting, i.e., the 
typical dual economic structure (rural versus urban sector) in industrializing countries. 
In section |2] the migration process is modelled within a statistical mechanics framework. 
Section E] presents the simulations and the main results. Finally, section El shows the 
concluding remarks. 



Let us consider an economic system formed by two sectors, one rural and the other urban. 
The main differences between these sectors are the sort of goods they produce, the tech- 
nologies used by firms and the framework of wage determination. Such a dual structure 
is typically used by the economic literature which investigates the rural-urban migration 
phenomena fT"3'] |H] . The basic features of the dual economy will be described in sub- 
section ITU and Subsection 12 .31 shows how the equilibrium macrostate of the economic 
system is determined. 

2.1 The urban sector 

The urban productive sector is formed by firms specialized in the production of manufac- 
turated goods. The output of the ith firm Ymi depends positively on both the amount of 
employed workers Nmi and the effort e, spent by each worker to perform his job. Based 
on the classical rural-urban migration theory we assume that the stock of capital 

during the analysis period is given. Supposing a standard geometrical functional form [7] , 
the production function of the manufacturing firm is the well-known Cobb-Douglas'^ 



where < a < 1 and Am > are parametric constants. 

By using the functional form originally suggested by Summers jH] and slightly modified 
by Romer |9| , the urban worker's effort can be defined as a function of the real wage paid 
by the manufacturing firm, the urban unemployment rate u and the alternative wage Wm, 
which is paid by other firms of the same sector. Then the effort function is given by 

^The modelling of employment and wage determination of the urban sector is based on the efficiency- 
wage approach. See chapter 10 of ref. and section V of ref. ^1] for further details. 
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where < rj < 1 and 6 > are parametric constants. 

The Zm manufacturing firms which form the manufacturing sector seek to maximize 
their real profits, measured in units of the manufacturated good, by choosing wages and 
employment freely. Given eq. and eq. (j2I) the real profit of ith manufacturing firm is 



Wmi - (1 - bu)Wr, 



^1 - hu)Wra J J 

The maximization condition of eq. (jHl) can be found using the first-order condition for 



a maximum, which result is 
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In equilibrium, all these firms choose the same wage jHlISl, i-e., Wmi = Wm = 
1, 2, 3, Zm)- Then, from equation (@)) the equilibrium urban unemployment rate is 

(6) 

By definition, the urban unemployment rate is the ratio between the number of unem- 
ployed workers and the urban population (A^„ — ZmNmi)/Nu, where Nu is the amount of 
workers localized in the urban sector. The previous definition must be consistent in each 
period to the equilibrium value of (jH)). The employment level of the manufacturing firm 
which obeys this consistency condition is obtained equaling the equilibrium in eq. (jHI) to 
the previous definition: 
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Taking eq. evaluated in the equilibrium, and eq. ^ and replace both in eq. ((H), 
the aggregated production of the manufacturing sector, Zml^mi, is given by 



V — N 
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(8) 



where a = A^^^-" [(^) (l - f ) 

By using eqs. (0), © and eq. (0 
facturing sector: 



one can obtain the equilibrium wage of the manu- 



^The second-order condition for a maximum is also satisfied. 
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Given the parametric constants that specify the technology and the effort sensitivity 
to wage, as well as the size of the manufacturing productive sector, it is possible to see 
that the equilibrium of the urban sector depends directly upon the urban population A^^. 



2.2 The rural sector 

In the rural sector the farm i produces an agricultural output Yai by employing an amount 
of workers Nai- The output is obtained by using a Cobb-Douglas production function [Zj 

Ya, = Aa{Na,)^, (10) 

where < < 1 and Aa > are parametric constants. We suppose that both the land 
endowment and the stock of capital of the farm are given during the period of analysis as 
assumed by refs. |3j and ^. 

Differently from the urban sector, farms are price-takers and the real wage is adjusted 
up to the point in which there is no unemployment in this sector This implies that 

the rural population will match the aggregated employment in the rural sector. Therefore, 
the equilibrium employment level of the farm i is 

iV„ ^ (U) 

where Za is the amount of farms which constitute the agricultural sector and N is the 
total population of the economic system. 

From eq. (fTU|) and eq. (fTT| the aggregated production of the rural sector, ZaYai, is 

= ^3 (iv - iv«)^ (12) 

where = ^a^i"*- 

Thus, the profit maximizing of the farms imply that the rural real wage expressed in 
units of the manufactured good becomes equal to the marginal product of agricultural 
labor in units of manufacturing good^: 

Wa = Up{N-N^)'f', (13) 

where ^4 = {Aa(j)/Z^-^) and p is the price of the agricultural good expressed in units of 
the manufactured good. 

Like in the urban sector, the equilibrium state of the rural sector depends on the urban 
population, as the size of total population of the economy is fixed. 

^This marginal product is the derivate of the production function, eq. (|10|l . with respect to Nai 
multiphed by p. 
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2.3 The macrostate of economic system 



As proposed by Harris and Todaro 13^, the terms of trade between the rural and urban 
sectors, measured by the price p, depend on the relative scarcity of agricultural and 
manufacturated goods. This can be measured by the ratio Ym/^a- The greater this 
ratio the greater will be the scarcity of agricultural good, which implies an increase of 
the agricultural good price in units of manufacturated good. Formally, given the urban 
population, the equilibrium relative price of the agricultural good is'' 



where p > and 7 > are parametric constants. 

Therefore, given the size of urban population, by using equations (jBEl) one can calcu- 
late the state of urban sector. Likewise, the rural sector state is determined by means of 
equations ()1HI14|) . The equilibrium state of both sectors will be modified if a migratory 
flux changes the population distribution of the economic system. 

3 Migratory process: a statistical mechanics approach 

As argued by Harris and Todaro j3]|3], individuals take their decisions of migrating or 
not by considering the differential of expected wages between their present sector and 
the sector they intend to go. However, other authors have taken into account additional 
reasons. Based on the formalization from statistical mechanics applied to socioeconomic 
phenomena [101 IH], in this section we propose an agent-based computational model to 
describe the rural-urban migratory process. This model is focused on the influence that 
individuals suffer in the reference group that they are included. The emergent properties 
will be analyzed taking into account the standard effect of labor allocation Harris- Todaro 
mechanism, which is based on the expected differential wages between sectors. This 
analysis will also be concerned on the effect of social neighborhood, often mentioned by 
other authors but not yet formalized. 

The main feature of the decision process is that each worker reviews his sectorial 
location after a period of time spent in that sector. We exclude, by assumption, the 
possibility that the worker may simultaneously supply his labor force to both sectors. 
Thus, only two choices are admitted: stay in the sector in which he was during previous 
periods or migrate. 

In order to model the migration process by allowing only discrete choices, each worker 
has its state defined by a-i G { — 1,+!}, where cxj = — 1 means that the worker is at the 
rural sector; otherwise, cXi = representing the urban sector. 

In our model, during the decision process, explicit and observable incentives are taken 
into account by each potential migrant. This is called a deterministic private utility 



^In the literature on rural-urban migration it is usual, because of analytical simplicity, to consider p 
constant 7_. This is true in the special case when 7 = in eq. (|14|) . 




(14) 



[TUj [TT| , given by 
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Ui = H{t)a,, (15) 

where H(t) = kue, > is a parametric constant and tUg is the expected urban-rural 
differential wage. The expected urban-rural wages in function of H{t) are specified as 
follows. 

Jobs are allocated at random when manufacturing firms are faced with more applicants 
than jobs available It nieans that in each time step all urban workers have the 

same probability to find an urban job. Under such a hypothesis, the term (1 — u) is the 
probability of an urban worker to obtain a job. Hence, (1 — u)Wm is the expected urban 
wage. Assuming that the rural wage is perfectly flexible there is no unemployment in 
the rural sector. Then, the probability to find a job in the rural sector is 1. Therefore, 
the rural wage Wa is the same as the expected wage in this sector. In sum, the expected 
differential of wage between urban and rural sectors is 

UJe = {I- u)Wm - Wa- (16) 

Besides, the worker i is also under the influence of other workers, his social neighbor- 
hood [12], denoted by nj. The measure of such influence, that is, the deterministic social 
utility [TU] [TT] , is given by 

Si = J^aiaj, (17) 

where J > is a parametric constant. The term J represents the interaction weight which 
relates the worker i's choice to the neighbor j's choice. This is assumed to be nonnegative, 
by representing the hypothesis that the worker seeks to conform to the behavior of his 
neighbors [llj. The interactions among neighbors are assessed in the workers' nearest 
neighbors or in the next nearest neighbors. 

Then, following references ^U] and ^1], we assume that payoff of worker i, which is 
his deterministic total (private and social) utility can be obtained replacing eq. (fT^ in 
eq. psp and summing with eq. (fT7|) : 

Hi = k[{l-u)Wm-Wa](Ti + J^(Tiaj. (18) 

Therefore, this system can be described by the well-known ferromagnetic Ising model, 
in the presence of an external time-dependent magnetic field: 

N 

n = -Hit) Y^cj^-JY. cjiaj. (19) 

i=l <ij> 

In each time step, each worker reviews his decision about the sectorial location with 
probability a, called activity [12]. Then, there is a part of the population that reviews 
their decisions and becomes potential migrants. 
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The potential migrant i becomes an actual migrant depending on the comparison 
between his deterministic total utility Tii and his non observable and idiosyncratic moti- 
vations fii, called random private utility JI]. The term fii represents the net difference 
between the random private utilities that the potential migrant assigns to the sector he 
intends to move and his present sector. 

In each period, if /ij > Tij, the potential migrant i becomes an actual migrant; other- 
wise, this does not happen. Supposing that /ij is a random variable logistically distributed 
^^[TT], the probability that the potential migrant effectively migrates is given by a cu- 
mulative distribution: 

where /3 > is a parametric constant that in this context measures the heterogeneity 
of workers concerning to the migration propensity. Equation ()20|) is a measure of the 
probability that a worker i, who is reviewing his location strategy, stays in the sector that 
he is localized at that time. The higher his deterministic total utility, eq. p8|l . the higher 
the probability that no change will take place. 



4 Simulation 

To carry out the simulation of the economic system described in the previous sections, 
each worker is placed in one of the sites of a square lattice. The state of each site (worker) 
is set as mentioned before: cXj = +1 for urban workers and Uj = —1 for rural ones. It 
is important to emphasize that the state of these sites represent the sectorial allocation 
of each worker, i.e., whether an individual is suppling his labor force in the urban or 
rural sector. It means that the coordinates of the lattice sites are not related to spatial 
distribution of workers. 

To set up the initial state of the system, all workers are randomly distributed in the 
lattice. At time t = there is the initial condition that 20% of the population is urban. 
In other words, initially, 20% of the sites will be assigned with cxj = +1 and the remaining 
80%, (Tj = —1. The reason for this initial distribution is because these are the values 
which have usually been observed in developing countries before the urbanization process 
initiates. 

The next step in the simulation is to calculate the equilibrium state variables of the 
urban sector, by using eqs. (jHlEl), and of the rural sector by using eqs. piH14p . Since 
the state variables of both sectors are known, it is necessary to define the amount of 
workers that will review their sectorial location, i.e. those one who will become potential 
migrants. To do this, it is assumed that the probability that a worker will become a 
potential migrant is given by the activity a, as defined by Stauffer and Penna |T^. All 
those selected as potential migrants will have their private utility calculated by eq. (fTSj) . 

In order to conclude the reviewing process, the probabihty defined in eq. ((20)) is 
assessed. Then, a random number rn e [0, 1] is generated from an uniform distribution. 
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Figure 1: Proportion of workers at urban sector as function of time for three different set 
of parameters J and k. Circles: ( J > 0, > 0); Squares: ( J = 0, A; > 0); Crosses: (J > 0, 
fc = 0). 

If rn > Pr, then the potential migrant becomes an actual migrant; otherwise, no change 
takes place. 

As soon as the potential migrants end their reviewing process, a new sectorial distri- 
bution is obtained. Knowledge of the new urban population allows the macrostate of the 
economic system to be reset. Therefore, the state variables of both sectors have to be 
calculated again. The whole procedure described above will be repeated as many times 
as we set in the simulation. The stopping criteria used by us is halting the simulation 
some steps after the moment when the system reaches equilibrium. 

Figure H shows the proportion of workers in the urban sector n„ = from now on 
called urban share, plotted in three different combination of the parameters J and k. It is 
necessary to remind that the parameters J and k adjust the instensity of the deterministic 
private utility, eq. (fT^ . and deterministic social utility, eq. (fTTj) . respectively. From top 
to bottom the set of parameters used in the plotting are (J > 0, /c > 0), (J = 0, /c > 0) 
and (J >0,k = 0). 

Firstly, consider the case (J = 0, A; > 0) plotted in Fig. [H In this case, the review 
conducted by the agents is guided only by the deterministic private utility, which in turn 
depends on the expected urban-rural difference of wages. As in models of classical theory 
of migration when the expected urban wage is higher than the rural wage, it implies 
in a continuous growth of the urban share, as well as a relatively fast convergence towards 
the equilibrium. 

Secondly, consider the case where both effects are taken into account, {J > 0, k > 0). 
Like the previous case, the rural-urban migratory process occurs again, however, the sys- 
tem reaches a higher value of the equilibrium urban share, though it takes more time for 
such outcome. This difference is caused by the parameter J > 0, what means that the 
influence of the social neighborhood is considered. To better understand this behavior, 
it should be reminded that the process of sectorial position revision depends on the de- 
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Figure 2: Rural share as function of per capita income in units of manufacturated good. 

terministic private utility and the social private utility. Then, when J > the influence 
of social neighborhood is being exerted, i.e., each worker attempts to adjust his choice 
according to the sectorial position of his neighbors. The existence of such an influence 
causes two different effects during the process of convergence towards equilibrium. In the 
first moment, when the neighborhood are mainly rural, the influence from neighbors slows 
the rural-urban migratory flux, increasing the time necessary to reach equilibrium. In the 
second moment, when the neighborhood become mainly urban, the influence reinforces 
the attraction from the high expected urban wage, leading to higher equilibrium urban 
share. 

Finally, we consider the case {J > 0, k = 0), with only neighborhood effects shown. In 
this case, the potential migrants consider only the sectorial position of the neighborhood 
and do not take into account the expected differential of wages. The pure effect due 
from neighborhood leads to the extinction of the urban sector. This is not an empirically 
important it has not been observed in developing economies. 

In Figure El another important feature caused by the migratory dynamics is the ex- 
pansion of per eapita income {Ym +pYa)/N. This result matches to the economic data in 
which in countries with high per capita income the proportion of the population living in 
rural area is low PP . 

In the initial state of the system the configuration was randomly set with 20% of the 
sites assigned cTj = +1, urban workers, and the rest cTj = —1, rural workers as shown in 
Fig. The final state of the dynamics by using (J > 0, A; > 0) can be visualized in 
Figure Eb- Now the infinite cluster is formed by sites cxj = +1 representing the urban 
concentration caused by the migratory process. Several others clusters are formed by sites 
CTi = -1. 

Figures 0] and El show the average magnetization m = 'Y^ai/N and the expected wages 
ratio Te = (1 — u)wm/u!a, respectively. Both figures are plotted as function of the ratio 
J/k {k kept constant) measuring the relative intensity between these parameters. 

Figure 0] has plotted in its vertical axis the average magnetization calculated during a 
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(a) 



Figure 3: Distribution of workers in the lattice. Black clusters represent urban workers 
(cjj = +1) and white clusters are rural workers (cTj = —1). Fig. (a) is the random 
distribution in the initial state of the system. Fig. (b) is the equilibrium distribution 
where clusters due the sectorial neighborhood can be seen. 

period after the system have reached equilibrium. To values of J / k less than the critical 
threshold the net magnetization is m = 0.4 representing an urban share about Uu = 0.70. 
By increasing the ratio J/ k after this critical threshold the system goes to a new regime, 
changing completely its net magnetization. 

Figure El is a plotting of expected wage ratio as function of J/k. To values J/k < 
9.0 the ration is Tf. = 1.0, what indicates that the expected urban wage and the rural 
wage converge to the same value. This property is known as Harris- Todaro equilibrium 
condition jSlCOIS]- Hence, in a economic system where internal migration occurs freely the 
absolute difference between the rural and urban wages can persist if workers consider the 
possibility of unemployment. After the threshold J/k > 9.0, Ve has its maximum value 
around 2.8 which shows that the urban expected wage is 2.8 times greater than the rural 
wage. Even having this ratio increasing the value of the worker private utility, eq. (jl5|) . 
the equilibrium of the system is m = —0.29, i.e., a rural concentration of 64.9%. The 
explanation of this outcome is that after a given threshold the values of J are in such a 
range that the social utility, eq. ()17|) . is many times higher than the private utility. In 
other words, in such range, it does not matter if the expected wage is attractive in the 
urban sector because the strongest factor in the migration decision is the influence of the 
neighborhood, i.e., agents tend to mimic the behavior of other agents. 

Simulations plotted in Figure 10 indicate that when the size of the lattice increases 
the equilibrium urban share Uu will change. For a given heterogeneity of the agents /3, 
there is a power law relating equilibrium urban share and the inverse of lattice size. This 
can be formalized in the expression below 
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Figure 4: Average magnetization m as function of ratio J/k. 
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Figure 5: Expected wages ratio as function of ratio J/k. 




where A and 6 are constants which have to be estimated. To carry out the estimation of 
these constants we evaluated a hnear regression of the log-hnear version of eq. (I2H)- In 
Table 1 one can find the estimation of the constants A and 6 based on data generated for 
five different values of (3. 

The estimation of the constants are approximately the same when using slightly dif- 
ferent values of agent heterogeneity (3. For example, by using any pair of constants A and 
6 from Table 1, the estimation of equilibrium urban share by eq. (|21|) is = 0.61 for an 
economy with 50 million of workers. 
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Figure 6: Log-log plot of equilibrium urban share n„ as function of the inverse of popu- 
lation size 
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Table 1. Estimates of parameters A and 9 for different values of (3. 

The effects of the ratio J/fc, together with size of population = L^, are shown in 
Fig. [7| The different values of equilibrium urban share are plotted in a grey scale. The 
first property observed in this figure is the existence of several phase states which depend 
on the values of J / k and A^. Each phase state is characterized by a constant equilibrium 
urban share. The topology of Figure [7| is in agreement with the results shown in Figs, 
m and El demonstrating that the properties of equilibrium macrostate depends on the 
combination of these parameters. 

In Figure IHl is plotted the equilibrium urban share as function of the parameter (3. For 
values of (3 tending to zero the equilibrium urban share tends to 0.5 (or m = 0), which 
implies in a null urban concentration (null average magnetization), even though there is an 
expected urban wage higher than the rural wage. In fact, eq. (j20p shows that the smaller (5 
the higher the idiosyncratic and non-observed proportion of the worker's behavior related 
to the migration propensity. If /5 = 0, the choices ai = +1 and cTj = — 1 have the same 
probability to occur being independent of the expected differential of wages. In sum, when 
the heterogeneity of the workers related to the decision of migration increases, the urban 
concentration will decline in the long run. On the other hand, when the heterogeneity 
of the agents decreases, i.e., /3 increases, the equilibrium urban share is invariable after a 
threshold. 
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Figure 7: Urban share rtu as function of the ratio J/k and the square lattice size L. 
Lighter areas correspond to higher n„ and darker areas to lower n^. 

The decision of migration is not taken simultaneously by all individuals. In order to 
simulate this behavior, the parameter called activity a is used. It gives the probability 
that a worker will review the decision about his sectorial location. More specifically, a 
represents the fraction of the population which will go through the reviewing process. 
This fraction of individuals is randomly selected and changes in each time step. In Fig. 
variation of a in different simulations shows that the time needed for the system to reach 
equilibrium is proportionally inverse to the value of the activity. Therefore, the time 
needed to reach the equihbrium state is strongly related to the amount of individuals 
which review their sectorial decision. 
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Figure 8: Equilibrium urban share n„ as function of the parameter (5. 
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Figure 9: Time to reach equihbrium as function of activity parameter a. 

5 Conclusion 

This paper has developed an agent-based computational model to analyse the rural-urban 
migration phenomena. The basic scenario was made of an economic system formed by two 
sectors, rural and urban, which differed in term of the goods produced, the production 
technology and the mechanism of wage determination. 

By assuming the sectorial migration decision as discrete choice in a miheu of de- 
centralized and non-coordinated decision-making, the rural-urban migration process was 
formalized as an Ising like model. The simulations showed aggregate regularities which 
indicates that decentralized migration decisions can lead to the emergence of equilibrium 
macrostates with features observed in developing economies. First, the simulation hav- 
ing an initial macrostate with population predominantly rural and expected urban wage 
higher than rural wage provoked a transitional rural-urban migratory dynamics, with con- 
tinuous growth of the urban share. This is a key feature of the phenomena called in ref. 
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[2] as urban transition. 

Second, simulations also showed that, during the rural-urban migration process, the 
reduction of the rural share takes place together with the increasing of per capita income 
of the economy. Such an inverse relation between rural share and per capita income is 
one of the most robust facts detected in economic statistics p. 

Third, the transitional rural-urban migratory dynamics converged towards an equilib- 
rium macrostate. The features of this transitional dynamics and equilibrium are sensitive 
to the relative weight between private and social effects (utilities) as well as the degree of 
heterogeneity of agents concerning the migration propensity. When the social interaction 
component is relatively stronger and below a critical threshold the transitional dynamics 
towards equilibrium is delayed and reaches a higher equilibrium urban share. With a high 
heterogeneity of agents, /? — 0, this generates the end of the pulhng force due the high 
expected urban wage what makes the system to reach an equilibrium macrostate with an 
urban share = 0.5. On the other hand, with a moderate heterogeneity of agents, /3 > 1, 
the equilibrium urban shares will be set in a empirically reasonable range (n„ > 0.6). 

The analysis shown in this paper suggests that a deeper investigation can still be 
carried out, which adopt alternative hypothesis mainly regarding the private and social 
utilities as well as other assumptions employed in our model. 
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